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Abstract In the typical spatial search problems solved by continuous-time 
quantum walk, changing the location of the marked vertices does not alter the 
search problem. In this paper, we consider search when this is no longer true. 
In particular, we analytically solve search on the “simplex of Km complete 
graphs” with all configurations of two marked vertices, two configurations of 
M + 1 marked vertices, and two configurations of 2(M + 1) marked vertices, 
showing that the location of the marked vertices can dramatically influence 
the required jumping rate of the quantum walk, such that using the wrong 
configuration’s value can cause the search to fail. This sensitivity to the jump¬ 
ing rate is an issue unique to continuous-time quantum walks that does not 
affect discrete-time ones. 

Keywords Grover’s algorithm • quantum search • spatial search • quantum 
random walk • multiple marked vertices 
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1 Introduction 


Schrodinger’s equation PEI is the fundamental equation of quantum mechan¬ 
ics that describes the evolution of a quantum state ij) in continuous time: 



= Hi!). 


We have set h = 1, and H is the Hamiltonian that characterizes the total 
energy of the system. For a particle of mass m in free space, the Hamiltonian 
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is simply the kinetic energy operator: 


H = 



where V 2 is the Laplace operator. 

The continuous-time quantum walk [3] is simply the discrete-space ana¬ 
logue of this. That is, the particle is confined to discrete positions in space, 
which can be expressed as the vertices of a graph, and to transitions expressed 
as the edges of the graph. Then the Laplace operator V 2 is replaced by its 
discrete version L = A — D, where A is the adjacency matrix of the graph 
(.Aij = 1 if two vertices i and j are adjacent, and 0 otherwise), and D is the 
diagonal degree matrix (Du = deg(?')). With this discrete substitution, the 
Hamiltonian for a continuous-time quantum walk is 


Hwalk "7 f- ? 

where we have also grouped the coefficients together into 7 , which is the jump¬ 
ing rate (i.e. amplitude per time) of the walk. Evolution by Sclirodinger’s 
equation with this Hamiltonian is precisely the definition of a continuous-time 
quantum walk US¬ 
As an algorithmic tool, continuous-time quantum walks have been used for 
a variety of applications, including evaluating boolean formulas [5], identify¬ 
ing graph isomorphism J 6 ], and performing universal computation [7]. In this 
paper, we focus on their application to search [3] on regular graphs. This leads 
to two changes to the Hamiltonian. First, since we assume that the graph is 
regular, the degree matrix D is simply a multiple of the identity matrix, so 
it can be dropped by rezeroing the energy or by factoring out a global, unob¬ 
servable phase. Thus we can use the adjacency matrix A instead of the graph 
Laplacian L. Second, we introduce a term that acts as a “Hamiltonian oracle” 
[Hj, which marks the k vertices to search for by potentials. With these two 
changes, the search Hamiltonian is 

H = - 1 A-Y j |i)<i|- (1) 

marked 

The goal is to evolve the system \ip(t)) by this Hamiltonian for as little time 
possible, from initially being in an equal superposition |s) over all the vertices, 

1 N 

\m) = \s) = -7= Ew. 

v i=l 

to a state that, when measured, collapses to a marked vertex with high proba¬ 
bility. Note that beginning in this state expresses our initial lack of knowledge 
of where the marked vertices might be, so it guesses each vertex with equal 
probability. Furthermore, this state is an eigenstate H w a ik, so the quantum 
walk alone does not change our information about where the marked vertices 
might be—it only changes when including the oracle, as in 0 - 
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Fig. 1 (a) The complete graph of N = 6 vertices, k = 2 of which are marked and denoted 
by double circles. Identically evolving vertices are identically colored and labeled, (b) The 
success probability as a function of time for search on the complete graph with N = 1024 
and k = 4. The black solid curve is when 7 = 7 C = 1/N, and the red dashed curve is when 
7 = 2 7 c = 2/N. 


As an example, let us consider the continuous-time quantum walk formula¬ 
tion of Grover’s algorithm mm ■ Since Grover’s algorithm solves the unstruc¬ 
tured search problem, one can move from any vertex to any other. So this is 
simply search on the complete graph of N vertices , an example of which is 


shown in Figure la 


By symmetry, the marked vertices evolve identically, as do the non-marked 
vertices. So we respectively group identically-evolving vertices together, as 
shown by identical colors and labels in Figure [Ta| 


» = ^E 

^ K ie red 


*). \ b ) = 


y/N — k 


E N>- 


iGblue 


Then the system evolves in a two-dimensional subspace spanned by {|a), 6)}, 
in which the search Hamiltonian Q is 

H = -J fc - 1 + i 

\s/k(N-k) N-k-1 J ' 

Following the analysis of mi, but generalized to k marked vertices, this has 
eigenstates 

. . , . , 1 - 7 IV ± AE 

“ |s> + 2VEV 

with gap in the corresponding eigenvalues Eq and E\ 

A E = E 1 -E 0 = V (1 - 7 ^) 2 + 4fcy. 


When 7 takes its critical value of j c = 1/N, the eigenstates are proportional 
to |s) ± \a) with an energy gap of AE = 2 y/k/N, so the system evolves 
from |s) to |a) in time n/AE = (n/2)^/N/k (see Section 3 of [12j for an 
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Fig. 2 Two equivalent ways to mark k = 2 vertices, denoted by double circles, on the 
complete graph of N = 6 vertices. 


explicit calculation). This is the quadratic speedup of Grover’s algorithm over 
a classical computer’s 0(N/k). As a check, Figure [Tb] shows the probability of 
measuring the quantum walker at a marked vertex as a function of time, and 
it reaches 1 at time (7r/2)-yi024/4 = 87 t « 25.133, as expected. 


From a straightforward calculation in Il3j (generalized to multiple marked 
vertices), 7 must be chosen within o(l/N 3 ' 2 ) of its critical value of 1/N for 
the algorithm to evolve from |s) to the marked vertices in time (tt/ 2)y/N/k for 
large N. This can be relaxed to evolve to the marked vertices with constant 
probability in time 0{^/N/k) if 7 is chosen within 0(l/iV 3 / 2 ) of its critical 
value of 1/N. If 7 is further away from its critical value than this, then the 
initial state |s) converges to an eigenstate of H for large N. So asymptotically, 
the system stays in |s) throughout its evolution, only picking up a global phase. 
This is shown in Figure lb with 7 = 27 c = 2 /N\ this is far enough from the 
critical value of 1/N such that the success probability plotted as a function of 
time converges to a flat, horizontal line for large N. 


Since Grover’s search problem is unstructured, the problem is unchanged 
no matter which k vertices are marked. For example, the two configurations 
in Figure [2] are equivalent (i.e., isomorphic). Since the location of the marked 
vertices does not change the structure of the problem, the algorithm (the 
critical jumping rate y c = 1/N, the evolution, the runtime of (tt/ 2)y/N/k, 
etc.) is unchanged. 

Even for search on non-complete graphs, which is called spatial search, it 
is possible to retain this property that the location of the marked vertices does 
not matter. The most common way is to search a vertex-transitive graph for a 
single marked vertex. This includes the hypercube [3], arbitrary-dimensional 
periodic square lattices [4], strongly regular graphs .14), the “simplex of com¬ 
plete graphs” na. and complete bipartite graphs |16| . For example, two pos¬ 
sible ways to mark k = 1 vertex on the simplex of complete graphs (defined 
more formally later) are shown in Figure |3j and they are clearly equivalent 
(i.e., isomorphic). Another way to make the arrangement of marked vertices 
irrelevant is by marking a cluster of vertices, such that moving the cluster 
leaves the search problem unaltered El- 
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(a) 


Fig. 3 Two equivalent ways to mark k 
simplex of complete graphs with M = 5. 



Even though the location of a unique (fc = 1) marked vertex on a vertex- 
transitive graph does not affect the search problem, when there are multiple 
marked vertices (k > 1), their locations generally do make a difference in 
spatial search. For example, there are five inequivalent ways to arrange k = 2 
marked vertices on the simplex of complete graphs, and they are shown in 
Figure [4j The simplex of complete graphs is the M-simplex with each of its 
M + 1 vertices replaced by a complete graph of M vertices, so it has a total 
of N = M(M + 1) vertices. In this paper, we explicitly solve spatial search 
by continuous-time quantum walk for these five configurations, plus four more 
with a greater number of marked vertices. In doing so, we show that different 
arrangements of marked vertices can dramatically change the critical jumping 
rate j c that is needed for the algorithm to succeed. 

This shows that the algorithm is dependent on the configuration of the 
marked vertices. Although such a dependence has been shown for discrete¬ 
time quantum walks | 18 l . this explicit demonstration seems new for continu¬ 
ous-time quantum walks. This result also highlights differences between the 
two approaches. While 7 C is primarily affected for continuous-time quantum 
walks, discrete-time quantum walks do not have this parameter. Instead, they 
are typically governed by a “coin” mm, and this coin is chosen in different 
ways to define a search problem m ■ For example, with the choice from 122| . 
one gets an algorithm that efficiently searches arbitrary-dimensional periodic 
square lattices with any configuration of two marked vertices m ■ With more 
marked vertices, the runtime can change 123] . even to the point of having no 
improvement over classically guessing j2¥j. With a different coin that yields 
precisely the phase flip in Grover’s algorithm mm, there are other excep¬ 
tional configurations that cause no improvement over classical as well [25]. 
Using Szegedy’s [26j method of defining a discrete-time quantum walk, one 
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Fig. 4 The five ways to distribute two marked vertices, indicated by double circles, on the 
simplex of complete graphs with M = 5. Identically evolving vertices are identically colored 
and labeled, and the labels indicate the subspace basis vectors that the vertices belong to. 
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can also obtain a quadratic improvement in search over a classical random 
walk’s “extended” hitting time when searching with multiple marked vertices 
m- This vast number of results for discrete-time quantum walks with multi¬ 
ple marked vertices dwarfs those for continuous-time quantum walks, of which 
this seems to be the first. 

We choose the simplex of complete graphs for our analysis because it has 
enough structure to reveal interesting properties, yet enough symmetry to 
be analytically tractable. It was first used in quantum search to prove that 
connectivity is a poor indicator of fast quantum search [lSj . and it has since 
been used to introduce a search algorithm that takes multiple walk-steps for 
each oracle query mi and to demonstrate faster search on a weighted graph 

m- 

In the next section, we solve spatial search on the simplex of complete 
graphs with two marked vertices, whose five possible configurations were shown 
in Figure [4j In doing so, we show that the first case’s critical jumping rate y c 
differs from the other four. Afterwards, we solve search with a greater number 
of marked vertices—two cases with M + 1 marked vertices and two cases with 
2 (M + 1) marked—showing that the critical jumping rate changes in a similar 
manner to the k = 2 case. This shows that for spatial search by continuous-time 
quantum walk, the critical jumping rate y c is dependent on the arrangement 
of the marked vertices. 


2 Two Marked Vertices 

The five possible configurations with k = 2 marked vertices were shown in Fig- 
ure[4j where the marked vertices are indicated by double circles. As shown, case 
(a) has both marked vertices in the same complete graph, and the remaining 
four cases (b), (c), (d), and (e) have them in different complete graphs. 

As with search with a single marked vertex [15], we get two-stage algo¬ 
rithms for each of these five cases, where the system evolves with one critical 
jumping rate 7 c i for some time, and then with a second critical jumping rate 
7 C 2 for a (likely different) amount of time. The detailed calculations for all five 
cases, including the generalization of the first case to any constant number of 
marked vertices in a single complete graph, are in Appendix [A] and the main 
results are summarized in Table [I] 

To give a sense of the calculations in Appendix[A]and explain the evolutions 
summarized in Table [l] consider the first case in Figure [4a] The system evolves 
in an 8 -dimensional subspace, independent of M, because there are only eight 
different kinds of vertices, as shown by the eight unique colors and labels in 
Figure |4a[ We group identically-evolving vertices into basis vectors |a), |5), 
..., \h) for the 8 D subspace. Then writing the search Hamiltonian ([I]) in 
this 8 D subspace, we find that for most values of 7 , the initial state |s) is 
asymptotically an eigenvector of H , which means the system does not evolve 
except for acquiring a global, unobservable phase. 
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Table 1 The five cases of search with 2 marked vertices, shown in Figure [ 4 ] with the 
subspace dimension and the two stage’s critical jumping rate ( 7 c ), runtime, and evolution. 


Case 

Dimension 


7c 

Runtime 

Evolution 

Figure 

4a 

8D 

3 

M 

+ 0 

( 1 > 
{.ms72) 

f M 3 / 2 

Is) -4 |6> 




1 

M 

+ 0 

f 1 > 
U 3 / 2 ) 


|6) — 4 |a) 

Figure 

4b 

4D 

2 

M 

+ 0 

( l \ 

5572^3/2 

Is) -»• 1 b) 




1 

M 

+ 0 

( 1 > 
U3/2 J 

f Vm 

|6) —4 |a) 

Figure 

4c 

8D 

2 

M 

+ 0 

( 1 \ 
Uf 5 / 2 ) 

W2^ 3/2 

Is) -+1&) 




1 

M 

+ 0 

( 1 > 
k m3/ 2 J 

f Vm 

\b) —■► l a ) 

Figure 

4d 

11D 

2 

M 

+ 0 

( 1 > 

^M 3/2 

Is) -»• 1 b) 




1 

M 

+ 0 

( 1 > 

| Vm 

|6) —4 |a) 

Figure 

4c 

13D 

2 

M 

1 

M 

+ 0 

+ 0 

( 1 > 
^ M 5 / 2 
( 1 ^ 
{m 3 / 2 J 

5F7 2M3/2 

f Vm 

\m) -> | b) + | h) 

|6) + | h) —> |a) + |d) 


1 

1 1 1 1 1 1 1 1 1 1 

T- 

1 

1 1 1 1 1 1 1 1 1 1 1 
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0.8 


- 
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- 
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- 
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0 
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0 
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Fig. 5 The success probability as a function of time for search on the simplex of complete 
graphs with M = 100 and k = 2 marked vertices, (a) Search for the first case using 
the correct 7 c ’s and runtimes for both stages of the algorithm (black solid) and using the 
incorrect values from the other cases (red dashed), (b) Search for the second case using 
the correct 7 c ’s and runtimes for both stages of the algorithm (black solid) and using the 
incorrect values from the first case (red dashed). 


For the system to evolve, we must choose the jumping rate 7 to be within 
o(l/M 5//2 ) of its critical value 7 cl = 3/M so that \s) experiences a “phase 
transition” [1], where it becomes supported by two eigenvectors of H instead 
of one. This can be found using degenerate perturbation theory nang, and 
it causes two eigenvectors of the Hamiltonian ([l]) to be asymptotically pro¬ 
portional to |g) ± | b) with an energy gap of A E = 6/M 3 / 2 . Since the equal 
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superposition state |s) is asymptotically | g) for large N (because the white 
g vertices in Figure [4a] overwhelming comprise most of the vertices in the 
graph for large N), near 7 = 7 C , the system evolves from |s) to | b) in time 


t\ = tt/AE = 7 tM 3 // 2 / 6 . From Figure 4a this means probability has now col¬ 


lected at the correct complete graph (ie., at the blue b vertices), but is not yet 
at the marked red a vertices within that complete graph. This is summarized 
in the first row of information in Table [T] 


For the second stage of the algorithm, we want to move the probability 
from |5) to the marked vertices |a). As shown in Appendix [A[ the Hamiltonian 
has two eigenstates that are asymptotically proportional to \b) ± a) when 
the jumping rate 7 takes the critical value 7 C 2 = 1/M with an energy gap 
of A E = 2 3 // 2 /-\/M. So the system now evolves from |5) to |a) in time £2 = 
tt/AE = 7 t\/M/ 2 3 ' /2 . This causes the system to evolve from the blue b vertices 
to the red a vertices, which are marked, accomplishing the search. This is 
summarized in the second row of information in Table [lj 


The success probability for the entire evolution of the algorithm is shown 
in Figure 5a as the solid black curve. For most of the time, the system is in 
the first stage of the evolution, where probability asymptotically builds up at 
| b). Any buildup in |a) during this time is due to higher-order contributions 
that are negligible for large N. Then in the second stage of the algorithm, 
the probability quickly shifts from | 6 ) to |a), indicated by the sudden spike in 
success probability in Figure [5a| 


Examining search for all five configurations in Table [T] reveals that the first 
configuration behaves differently from the other four—it has a different critical 
jumping rate for the first stage of the algorithm 7 c i, and the runtime of both 
stages is different. Of these, the important difference is the critical jumping 
rate, which truly is different because 3/M and 2/M are separated by a gap of 
more than o(l/M 5 / 2 ). Thus evolving by the wrong configuration’s value would 
cause the system to asymptotically stay in its initial state |s) throughout 
its evolution, only acquiring a global, unobservable phase. This is shown in 
Figure [5a[ where the dashed red curve is search for the first configuration, but 
incorrectly using the critical 7 ’s and runtimes from the other cases. Note that 
the probability stays small, and any buildup from higher-order corrections 
vanishes for large N. Similarly, Figure [5b] shows search with vertices in the 
second configuration in Figure [4b] -with the correct critical 7 ’s and runtimes, 
success probability builds up as desired, but with the wrong values from the 
first configuration, it fails to build up. Note that the second, third, fourth, and 
fifth configurations corresponding to Figures [4bJ[4cJ[4dJ and[4e[ have identical 
critical jumping rates and runtimes. So while rearranging marked vertices can 
change the critical jumping rate substantially, it could also do nothing. 


Given the extensive calculations in Appendix [A] that are needed to de¬ 
rive the results in Table [TJ manually working through all configurations of k 
marked vertices on the simplex of complete graphs is impractical using the 
current method. We do analyze in the next section, however, four additional 
configurations with a greater number of marked vertices, two of which show 
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Fig. 6 Two ways to distribute M + 1 marked vertices, indicated double circles, on the 
simplex of complete graphs with M = 5. Identically evolving vertices are identically colored 
and labeled, and the labels indicate the subspace basis vectors that the vertices belong to. 




(a) (b) 

Fig. 7 Two ways to distribute 2(M + 1) marked vertices, indicated double circles, on the 
simplex of complete graphs with M = 5. Identically evolving vertices are identically colored 
and labeled, and the labels indicate the subspace basis vectors that the vertices belong to. 


that moving marked vertices again changes the critical jumping rate, and two 
that show no change. 


3 Larger Examples 

Let us consider four examples with a larger number of marked vertices. The 
first two are shown in Figure |6j and they are just two of the many ways to 
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Table 2 Some cases of search with multiple marked vertices with the subspace dimension, 
critical jumping rate (7c), runtime, and evolution. 


Case 

Dimension 

7c 

Runtime 

Evolution 

Figure 

Figure 

6 a 

6b 

3D 

7D 

m + °(j^m) 

1+ M +°( M 3/ 2 ) 

fVM 

fVM 

|6) —> \a) 

1 9) -> |6> 

Figure 

Figure 

7a 

7b 

2D 

3D 

M + 0 (jj372j 

M + 0 ( M 3/2) 

575^ 

275^ 

|6) —> | a) 

|fe) —> |a) 


arrange M +1 marked vertices on the simplex of complete graphs; in subfigure 
(a), they are distributed so that each complete graph has one marked vertex, 
and in subfigure (b), all the marked vertices are clustered in a single complete 
graph, except for one. The other two configurations are shown in Figure[7j and 
they are just two of the many ways to distribute 2 (M + 1) marked vertices so 
that each complete graph has two marked vertices. 

Unlike search with a single marked vertex m or two marked vertices in 
the previous section, search on these graphs are single-stage algorithms. The 
detailed proofs are given in Appendix [Bj and they use the same techniques as 
the configurations with two marked vertices. The results are summarized in 
Table [2} and they show that the critical jumping rate q c for the two configu¬ 
rations in Figure [b] differ substantially enough that the search algorithm will 
fail (i.e., the system will asymptotically stay in its initial state) if the wrong 
7 C is used. In Figure [7J however, the critical jumping rate y c and runtime are 
the same. So again, we see that rearranging the marked vertices can change 
the critical jumping rate substantially, but it also might not change it at all. 


4 Conclusion 

We have shown that when there are multiple marked vertices, their configura¬ 
tion on a graph can affect the critical jumping rate y c of the continuous-time 
quantum walk. Previous work on search by continuous-time quantum walk has 
avoided this effect by restricting search to a single marked vertex on vertex- 
transitive graphs, or by clustering marked vertices together. This highlights a 
difference between continuous- and discrete-time quantum walks as they are 
used for search, since discrete-time quantum walks do not have the jumping 
rate as a parameter. 

This work leaves open how general configurations of multiple marked ver¬ 
tices affect search on the simplex of complete graphs, since our analysis only 
examined nine specific arrangements. As a curious observation, our results fol¬ 
low a pattern: the critical jumping rate 7 C is affected by the number of marked 
vertices in each complete graph, not by their arrangement within them. That 
is, going from Figure to Figure [6j:> moves marked vertices across com- 
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plete graphs, changing the number within each, which changes 7 c . On the 
other hand, going from Figure [TJi to Figure [?]) moves marked vertices within 
complete graphs, and this rearrangement within complete graphs makes no 
difference to 7 C . This behavior is consistent with two marked vertices in Fig¬ 
ure^ and it might stem from each complete graph being sufficiently connected 
within itself that the arrangement of marked vertices within them does not 
matter. Whether this pattern holds in general is a subject of further inves¬ 
tigation, as is how continuous-time quantum walks search other graphs with 
multiple marked vertices. 

Acknowledgements Thanks to Andris Ambainis for useful discussions. This work was 
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A Details for Two Marked Vertices 

In this appendix, we employ degenerate perturbation theory ECUS to find the critical 
7 ’s and runtimes for search with two marked vertices, of which there are five cases, as 
summarized in Table [l] 


A.l Two Marked, Case 1, Generalized to Constant Marked Vertices 

Instead of having just 2 marked vertices in a single complete graph, we generalize the problem 
to k constant marked vertices. Even with this generalization, the system still evolves in an 
8 -dimensional subspace, as shown in Figure |4a| spanned by 





X 10, 


i G magenta 
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Fig. 8 Apart from a factor of — 7 , (a) the Hamiltonian for the first case of search on the 
simplex of complete graphs with k = 2 marked vertices, (b) the leading-order terms for 
the first stage of the algorithm, and (c) the leading-order terms for the second stage of the 
algorithm. 


In this subspace, the search Hamiltonian 0 is 


H = — 



0 7 M k (k - 1 ) 0 



where = M — k and = M — k — 1. 

Using the diagrammatic approach in m as a guide, this Hamiltonian can be visualized 
as a graph with eight vertices, as shown in Figure [ 8 ^.. For the first stage of the algorithm, 
the leading-order Hamiltonian i/(°) can be visualized as shown in Figure [ 8 |d, where we have 
excluded edges that scale less than y/~M. From this, the eight eigenvectors of are easily 
seen: two are linear combinations of | a) and | 6 ), three are linear combinations of |c), | d), and 
| h), and the final three are linear combinations of |e), |/), and | g). They correspond to the 
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eigenvectors of 


jr (0 > = - 7 ( 

ab 1 fkM M 


H 


(0) _ 
cdh 


( 0 


Vm o 

= -7 I Vm m V M ( k - !) 

\ 0 y/M{k - 1 ) 0 

o o Vm 

0 0 VkM 

Vm VkM m 


Jf<°> = -7 
efg ' 


Since |s) « |g), and we want probability to move towards the marked vertices |a), we want 
to choose 7 so that a linear combination of |e),|/), and | g) is degenerate with a linear 
combination of |a) and | b). In particular, the eigenstates that we want to be degenerate are 


7>- 


2 VkM 


Vm + V4 + 4 k + M Vm (Vm + V4 + 4k + m\ 

with corresponding eigenvalue 

E U = V1 (m + VmV 4 + 4fc + m) 

and 


l/> + Is) 


1 — M 7 + y/l — 2 M 7 + 4 /CM 7 2 + M 2 7 2 
2\/ r kMj 


I a) + |6> 


with corresponding eigenvalue 

E v =: i (-1 -My- Vl - 2 M 7 + 4k My 2 + M 2 7 2 ) . 

Written this way, u and v are unnormalized, whereas \u) and |u) are their normalized 
versions. These eigenstates are degenerate when 7 takes its critical value of 


7ci = 


-M + VMV 4 + 4 k + M 1 + fc 


2 M 


M 


The perturbation H^\ which restores terms of constant weight, causes certain linear com- 
binations 

ol u \u) + a v \v) 

of these states to be eigenstates of + H I 1 ' |2lll4|. The coefficients a u and a v can be 
found by solving 

/Huu H uv \ /O-uA _ £ {Ot-x 

\H VU H vv J yOi-u J 

where H U v = (u\H^ + H^\v), etc. Solving this, the perturbed eigenvectors for large N 
with their corresponding eigenvalues are 

1 „ v „ . k 2 +2k + 1 k + 1 

— (\u) + \v)), E--(k + 1)+ — jjpya 

1 „ , , . , , k 2 + 2k + 1 k + 1 

-(\u)-\v)), E--{k + 1)+ — + J ^ 75 

Since |u) « \g) and |u) ~ |6) for large N, the system evolves from |s) ~ | g) to |6) in time 
n/AE , which is 

tM 3 / 2 


1 1 = 


2(fc + 1) 
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Diagrammatically, the perturbation H restores edges of constant weight in Figure]^, and 
probability flows between | g) and | 6 ) since they are the most dominant terms. 

Using the approach of Section VI of m, if 7 is within e of its critical value of 7 c i £3 
(1 + k)/M , then the eigenvalues of |u) and |u) now include leading-order (in e) terms —eM. 
In the perturbative calculation, this introduces terms scaling as eM due to H uu and H vv , 
so for this to not influence the energy gap @(1/M 3 / 2 ), we must have eM = o(l/M 3 / 2 ), or 
e = o(l/M 5 / 2 ). Thus for the first stage of the algorithm to asymptotically evolve from |s) 
to | 6 ), we require 7 = 7 c i + o(l/M 5 / 2 ). Note if we relax this to evolve to | b) with constant 
probability, then 7 = 7 c i + 0(1/M 5 / 2 ) suffices. 

For the second stage of the algorithm, we take the leading-order Hamiltonian if (°) to 
only include edges of weight <9(M), and its diagram is shown in Figure [sj:. From this, the 
eight eigenvectors of if (°) are simply the basis vectors |a), | 6 ), ..., | h) with corresponding 
eigenvalues —1, — 7 M, ..., 0. When 7 takes its critical value of 


the eigenstates | a), |6), |d), and | g) of ifare degenerate with eigenvalue —1. Then the 
perturbation H^\ which restores terms <9(\/M), causes certain linear combinations 


a a \a) + OL b \ b) + oL d \ d) + a g \g) 


of these states to be eigenstates of if+ H^\ The coefficients a a , ot bi a d> and a g can be 
found by solving 


' EL aai E a b E ad H a g' 


( <0 



Eba E bb H-bd Ebg 



= E 


Eda Hdb Hdd Edg 





\Hga Hgb Hgd Egg) 


\ a sJ 


\ a sJ 


where H ab = (a|if(°) + H^\b), etc. Solving this, the perturbed eigenvectors with their 
corresponding eigenvalues are 


IV>o} 


|V>3> 


-^(|6)-(-|a)), E 0 - -1 

iVh) = \d), E 1 = -l 
l^ 2 > = \g), E 2 = -1 

-4(|b>-|a», Ea = - 1 




So the system evolves from |6) to |a) in time tt/AE: 


7 r M 

t2= 2VT- 


Diagrammatically, the perturbation if^ 1 ) restores edges (9(\/M) in Figure [SJ:, and probabil¬ 
ity flows between | 6 ) and |a). 

We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 C 2 = 1/M—a straightforward calculation shows that it must be within o(l/M 3 / 2 ). 
For k = 2 marked vertices, the critical 7 ’s and runtimes are 


7d = 


M' 


1 1 — 


T M 3 / 2 
_ 6 ’ 


7c2 = 


M’ 


t2 = 


M 

2~’ 


all of which are in agreement with Table [T] 
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Fig. 9 Apart from a factor of —7, (a) the Hamiltonian for the second case of search on 
the simplex of complete graphs with k = 2 marked vertices, (b) the leading-order terms for 
the first stage of the algorithm, and (c) the leading-order terms for the second stage of the 
algorithm. 


A.2 Two Marked, Case 2 


As shown in Figure [ 4 b| the system evolves in a 4 -dimensional subspace spanned by 

1 


‘ 7! £ '*>• 


It). 


iEred 

i 


V 2 ( M ~ 1) 


|6> V2(M-1) ie £ n |l> ’ 

l9>= ^/{m -i){W^T) ^J l) ’ 


where the labels have been chosen this way to match the behavior of the vertices in the first 
case in Figure [ 4 a] In this subspace, the search Hamiltonian § is 


H = —7 


/ l+A 0 0 \ 

y/M - 1 M — 2 1 0 

0 1 1 y/2 (M - 2) _ 

\ 0 0 -J 2 {M- 2) M -2 J 


This Hamiltonian can be visualized as shown in Figure [ 9 jt. For the first stage of the 
algorithm, the leading-order Hamiltonian [f uy> excludes edges that scale less than \fM , and 
it can be visualized as shown in Figure]!)]]. The two eigenstates of /A 0 ) that we want to be 
degenerate are 

VM-yfM+~8 
U =-7=-|e) + Is) 


2 V 2 


with corresponding eigenvalue 


and 


E u = [M + yfMy/M + 8^ 

— 1 + M7 — yjl — 2M7 + 4 M 7 2 + M 2 7 2 


I a) + |6> 


with corresponding eigenvalue 

E v = i (-1 - M'y- y/l - 2M7 + 4 My 2 + M 2 7 2 ) . 
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These are degenerate when 7 takes its critical value of 

= -M + VWM + 8 = 2__4 

2M M M 2 K J 

The perturbation H^\ which restores terms of constant weight, causes certain linear com¬ 
binations a u \u) +a v \v) to be eigenstates of -|-Doing the perturbative calculation 

to find the coefficients (as in the first case), the perturbed eigenstates for large N are 

T(|«) + |„)), B— 

75<l“> 

Since | u) « | g) and |u) « \b) for large N, the system evolves from |s) « | g) to | 6 ) in time 
n/AE : 

tM 3 / 2 


1 1 = 


4y/2 


Diagrammatically, the perturbation H ^ restores edges of constant weight in Figure]^), and 
probability flows between |g) and | 6 ) since they are the most dominant terms. 

As in the first case in the precious section, we can use the method of m to find how 
precisely 7 must be chosen to its critical value 7 c i = 2/M —a straightforward calculation 
shows that it must be within o(l/M 5 / 2 ). 

For the second stage of the algorithm, we take the leading-order Hamiltonian to 
only include edges of weight (9( M ), and its diagram is shown in Figure When 7 takes 
its critical value of 

7c2 = i’ 

the eigenstates |a), | 6 ), and | g) of H are triply degenerate. Then the perturbation 
which restores terms <9(\/M), causes certain linear combinations a a \ci)-\-&b\b)-\-Oig\g) of them 
to be eigenstates of H Doing the perturbative calculation to find the coefficients 
(as in the first case), the perturbed eigenstates for large N are 

+ e = - 1 ~7m 

\g), e = -i 

^ (|6) - |a)) ’ e = - 1 + 7m ■ 


So the system evolves from | 6 ) to |a) in time n/AE: 


7 ry/M 


Diagrammatically, the perturbation restores edges &(y/M) in Figure^, and probabil¬ 
ity flows between | 6 ) and |a). 

We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 c 2 — 1 /M —a straightforward calculation shows that it must be within o(l/M 3 / 2 ). 
These 7 c ’s and runtimes are in agreement with Table [I] 
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Fig. 10 Apart from a factor of — 7 , (a) the Hamiltonian for the third case of search on the 
simplex of complete graphs with k = 2 marked vertices, and (b) the leading-order terms for 
the first stage of the algorithm. 


A.3 Two Marked, Case 3 


As shown in Figure [4c| the system evolves in an 8 -dimensional subspace spanned by 

1 


l“> = 4 E !*>’ 

VZ ie red 

|6> = 1 y i»), 

V2(M-2) ^ 

| C ) = ^, E W. 


\d) = 


iGyellow 

i 


yJM - 2 


E i*>' 


|e> = 

l/> = 

Iff) = 


E i®> 


s/M - 2 


iG brown 

i 


._= y i i), 


my e !*)• 


iG magenta 


iGgray 


Note there are no h type vertices. Instead, there’s a new type, which we call i. In this 
subspace, the search Hamiltonian 0 is 

/ i \/M2 1 0 0 0 0 1 \ 

V m 2 A/3 00100 v m 2 

0 
0 
0 
0 
0 

1 / 

This Hamiltonian can be visualized as shown in Figure [To},. For the first stage of the 
algorithm, the leading-order Hamiltonian H (°) excludes edges that scale less than y/~M, 
and it can be visualized as shown in Figure 01 - As with the last two cases, there are two 
eigenvectors that we want to be degenerate. The first is 

2 \/2 2 

y/M + y/M + 12^ + y/M + y/M + 12^ + ^ 
with corresponding eigenvalue 


H = - 7 


1 
0 
0 
0 
0 

V 1 


0 

0 

1 

0 

0 

VM2 


1 V2M2 
/ 2 M 2 M 3 
0 0 


where M 2 = M — 2 and M 3 = M — 3. 



/ 2 M 3 >/ M 3 M 3 
0 0 0 




= E (^ + y/M(M + 12)) . 
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The second eigenvector is messy, but can be approximated nicely. The leading-order Hamil¬ 
tonian corresponding to |a), |fo), and \i) is 


H 


(0) = - 7 
a.b.i 1 


( i Vm o \ 

Vm m Vm 

y 0 VM 0 J 


The eigenvalues A of this satisfy the characteristic equation 

-A 3 - ( 7 M + 1)A 2 + 7 M(2 7 - 1)A + 7 2 M = 0. 


When 7 takes its critical value of 


2 _ 6 _ _36_ 

~M~ ~M 2 + M 3 ’ 


one of these eigenvalues, which we will call E v , and E u both equal — 2 — 270/M 3 + 0(1/M 4 ), 
making them approximately degenerate. To find the corresponding eigenvector v, we use the 
first and third lines of the eigenvalue equation so = E v v: 


-7 


( \ Vm o \ 
Vm M s/M 
y 0 s/M 0 ) 



This yields 


-ls/M 

v = - a) • 

E v + 1 1 ' 




The perturbation which restores terms of constant weight, causes certain linear com¬ 

binations a u \u) -\-a v \v) to be eigenstates of H -| -H^\ Doing the perturbative calculation 
to find the coefficients, the perturbed eigenstates for large N are 




E = - 2 - 


2s/2 

M 3 / 2 


E = -2 + 


2s/2 

M 3 72' 


Since |n) Ri |g) and \v) 
tt/AE: 


16) for large N, the system evolves from s) 


11 


n rM 3 / 2 
4v/2 


| g) to | b) in time 


Diagrammatically, the perturbation restores edges of constant weight in Figure fTTT] ;. 

and probability flows between | g) and \b) since they are the most dominant terms. 

We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 c 1 = 2/M —a straightforward calculation shows that it must be within o(l/M 5 / 2 ). 

The second stage of the algorithm is similar to the previous two cases, where we take 
the leading-order Hamiltonian to only include edges of weight &(M). When 7 takes 
its critical value of 

1 


7c2 = 


Af* 


then |a) and \b) are degenerate eigenvectors (among others) of The perturbation H^ 

restores terms of order 0(vM), which causes probability to flow between | 6 ) and | a). Doing 
the calculation, we find eigenstates of the perturbed system that are proportional to | 6 ) zb \a) 
with eigenvalues — 1 =p l/y/~M , so the system evolves from \b) to |a) in time tt/AE: 


7 ry/M 


As before, a straightforward calculation using the method of m shows that 7 must be 
chosen within o(l/M 3 / 2 ) of its critical value 7 C 2 = 1/M. 

These 7 c ’s and runtimes are in agreement with Table [I] 
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A.4 Two Marked, Case 4 


As shown in Figure [4d| the system evolves in an 11-dimensional subspace spanned by 


|a) = 

b) - 

|c> = 
\d) = 
\e) = 

I/> = 


1 


E 10. 

iEred 
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l<7> = 
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V (M — 3)(Af — 4) ie £ it J l> ’ 


__ Y" |i), 

V2(M-3) ie tr ue 




li), 
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E io. 

iEyiellow 
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__ y i i), 

V2(M-3) iem ^ enta 
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10 = 


li> = 


V2 


t/2 


iE orange 


E io. 

iEgray 

E i*>> 

iEpink 


'*> - 73 £ '«■ 


iGlime 


In this subspace, the search Hamiltonian § is 


H = -7 



1 1 

V M 3 V-^3 

0 
0 
0 
0 
0 
0 
1 
0 
1 


0 ' 

0 

1 

s/M/, 

0 

0 

0 

1 

0 

1 

0 , 


where M 3 = M — 3 and M 4 = M — 4. 

This Hamiltonian can be visualized as shown in Figure m- For the first stage of the 
algorithm, the leading-order Hamiltonian H excludes edges that scale less than \/M, 
and it can be visualized as shown in Figure E3>- As with the last two cases, there are two 
eigenvectors that we want to be degenerate. The first is 


2v/2 


s/M + /W- 

with corresponding eigenvalue 

E,. = 


16 


l e ) ' 


2v/2 


s/M + s/M + 16 


l/> + lfl>. 


-7 


M + s/Ms/M + 16) . 


The second eigenvector is messy, but can be approximated nicely. The leading-order Hamil¬ 
tonian corresponding to |a), | b), |i), and | j) is 


H (0 1 ■ ■ = -7 

a,b,z,j 1 


( i s/M 0 0 \ 

s/M M s/M s/M 

0 s/M 0 0 

\ 0 s/M 0 0 ) 

The eigenvalues A of this satisfy the characteristic equation 

A 4 + (1 + M 7 )A 3 + M/l - 37 )A 2 - 2M 7 2 A = 0. 
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M 


Fig. 11 Apart from a factor of — 7 , (a) the Hamiltonian for the fourth case of search on 
the simplex of complete graphs with k = 2 marked vertices, and (b) the leading-order terms 
for the first stage of the algorithm. 


When 7 takes its critical value of 


2 _ 8 64 

7cl “ M ~~ M 2 + M3’ 

one of these eigenvalues, which we will call E v , and E u both equal —2 —640/M 3 — 0(1/M 4 ), 
making them approximately degenerate. To find the corresponding eigenvector v, we use the 
eigenvalue equation i .v = E v v: 


-7 
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0 

rH|f 




V 

s/M M s/M s/M 


Vb 

= E v 

Vb 

O 

<• 

s: 
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0 


Vi 


Vi 

\ 0 s/M 0 0 J 


\ v i) 


\ V J 


This yields 


v = 


—7 y/~M 
1 + E 


I °) + | b) 



7 V M 
E 


I j). 


The perturbation H^\ which restores terms of constant weight, causes certain linear com¬ 
binations a u \u) +ol v \v) to be eigenstates of H -| -H^. Doing the perturbative calculation 
to find the coefficients, the perturbed eigenstates for large N are 


-4(i u )-w))- 


E — —2 — 


2V2 

M372 


E = - 2 + 


2V2 

M 3 /2 ' 


Since \u) Ri r/} and ?;} « | 6 ) for large N, the system evolves from s) ps | g) to \h) in time 
■rr/AE: 

ttM 3 / 2 

Diagrammatically, the perturbation restores edges of constant weight in Figure ink 
and probability flows between |^) and | b) since they are the most dominant terms. 
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We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 d = 2/M—a straightforward calculation shows that it must be within o(l/M 5 / 2 ). 

The second stage of the algorithm is similar to the previous three cases, where we take 
the leading-order Hamiltonian H to only include edges of weight <9(M). When 7 takes 
its critical value of 


then |a) and | 6 ) are degenerate eigenvectors (among others) of The perturbation H^' 

restores terms of order <9(\/M), which causes probability to flow between \b) and |a). Doing 
the calculation, we find eigenstates of the perturbed system that are proportional to \b) =b |a) 
with eigenvalues —1 =p 1/vM, so the system evolves from | b) to |a) in time n/AE: 


7 ry/M 


As before, a straightforward calculation using the method of m shows that 7 must be 
chosen within o(l/M 3 / 2 ) of its critical value 7 C 2 = 1/M. 

These 7 c ’s and runtimes are in agreement with Table ^ 


A.5 Two Marked, Case 5 


As shown in Figure [4e| the system evolves in a 13-dimensional subspace spanned by 
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In this subspace, the search Hamiltonian 0 is 
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Fig. 12 Apart from a factor of — 7 , (a) the Hamiltonian for the fifth case of search on the 
simplex of complete graphs with k = 2 marked vertices, and (b) the leading-order terms for 
the first stage of the algorithm. 


where M 2 = M — 2 and M 3 = M — 3. 

This Hamiltonian can be visualized as shown in Figure ESt. For the first stage of the 
algorithm, the leading-order Hamiltonian H excludes edges that scale less than y/M, and 
it can be visualized as shown in Figure EU 3 - The initial equal superposition state |s) is 
approximately | m) for large N , and we want it to evolve to the marked vertices |a) and | d). 
So we will need leading-order eigenstates that are approximately each of these to be triply 
degenerate. The first is 

(|e> + |fc> + | 0 ) + | m) 


\fM + sj 12 + M 
with corresponding eigenvalue 


E u = —^ (m + sjM{M + 12)) . 


Note this is t he s ame eigenvalue as E u from Case 3. For the other two leading-order eigen¬ 


states, Figure |12 Jd reveals that i and h are identical with a ~ d, b ~ h, and i ~ c, 
so their correspond ing eigenstates are always degenerate. Furthermore, they are identical to 

(o) - rtl 


( 0 ) 


H 


a,b,i 


from Figure 1 Op from Case 3. So the eigenvectors and eigenvalues carry over: 




—jy/M -Jy/M 

l“> + l b > + —^—W, 


E v + 1 
—7 \/M 
Eyj + 1 


\d) + \h)+^^\c), 


including the critical 7 


6 


36 


7cl M M 2 + M 3 ’ 


at which the eigenvalues E u , E v , and E w all equal —2 — 270/M 3 + 0(1/M 4 ), making them 
approximately degenerate. 

With the perturbation which restores terms of constant weight, we have the same 

behavior and runtime 

7 rM 3 / 2 
tl = - -j=- 

4\/2 
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as Case 3, except the system evolves from |s) « | m) to | 6 ) +1 h). So the probability gets split 
between the two paths. Diagrammatically, the perturbation H B) restores edges of constant 
weight in Figure [T2fc , and probability flows from | m) to \b) and | h) since they are the most 
dominant terms. 

We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 C 2 = 1/M—a straightforward calculation shows that it must be within o(l/M 5 / 2 ). 

The second stage of the algorithm is similar to the previous cases, where we take the 
leading-order Hamiltonian H to only include edges of weight <9(M). When 7 takes its 
critical value of 


then |a) and | b) are degenerate eigenvectors, as are | d) and | h), of . The perturbation 
H restores terms of order <9(\/M), which causes probability to flow from | 6 ) to |a) and 
from | h) to | d). The runtime from Case 3 carries over: 


7r\/M 


As before, a straightforward calculation using the method of m shows that 7 must be 
chosen within o(l/M 3 / 2 ) of its critical value 7 C 2 = 1/M. 

These 7 c ’s and runtimes are in agreement with Table [I] 


B Details for Larger Examples 

In this appendix, we employ degenerate perturbation theory GKOES] to find the critical 
7 ’s and runtimes for search with a larger number of marked vertices, the results which are 
summarized in Table |2] 


B.l One Marked Per Complete 

As shown in Figure | 6 ^, the system evolves in a 3-dimensional subspace spanned by 


|a) = 7 1 Y |i) 
s/M + 1 1 ' 

1 




| 6 > = , ' ^ Y |*>, 

|c> = , 1 V |i>. 

v/M + 1 


y/M 

In this subspace, the search Hamiltonian 0 is 


£E yellow 


1 yJM- 2 2 

H = -7 ( CM-2 M - 2 y/M- 2 ) . 

2 y/M -2 0 


This Hamiltonian can be visualized as shown in Figure [13^ . The leading-order Hamil¬ 
tonian H excludes edges that scale less than M, and it can be visualized as shown in 
Figure [Hfo . Clearly, the eigenstates of this are |a), 16), and |c) with corresponding eigen¬ 
values — 1, — 7 M, and 0. Since |s) ~ | 6 ), we choose 7 so that | 6 ) and |a) are degenerate, 
i.e., 















26 


T. G. Wong 


M -2 7 M 



Fig. 13 Apart from a factor of —7, (a) the Hamiltonian for the first case of search on the 
simplex of complete graphs with k = M + 1 marked vertices, and (b) the leading-order 
terms. 


The perturbation H^\ which restores edges of weight causes certain linear com¬ 

binations <y. a \o) + ot b \b) to be eigenstates of H + H (■*■). The coefficients can be found by 
solving 

fHaa H ab \ f CfccA _ rp f&a\ 

\H ba H bb ) \a b ) -*{<*)' 

where H ab = (a\H^ + H^\b), etc. With 7 = 7 C , this yields eigenstates and eigenvalues 

\ipo) = ^g(l 6 ) + l°»> E 0 

\1>i) -^06) - |o)), E 1 

So the system evolves from |s) « \b) to |a) in time n/AE: 

7 T\/ M 
U = ——• 

2 

Using the method Section VI of m , if 7 is within e of its critical value of 7 C = 1 /M, 
then the eigenvalue of \h) is now —7 M = — 1 — eM. In the perturbative calculation, this 
introduces a leading-order (in e) term eM due to Ih>h- For this to not influence the energy 
gap 6>(l/\/M), we require eM = o(l/\/M), or e = o(l/M 3 / 2 ). Thus for the algorithm to 
asymptotically evolve from \b) to |a), we require 7 = 7c + o(l/M 3 / 2 ). 

This 7 C and runtime are in agreement with Table | 2 | 
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B.2 Fully Marked Complete, Plus One 


As shown in 


Figure [6)3, the system evolves in a 7 -dimensional subspace spanned by 


|a) = | red), 


|6> 


1 

sJM - 1 
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iGblue 


|c> = I yellow), 


M> 
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y/M - 1 
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l e > 

l/> 

la> 


1 

y/M - 1 
1 

y/M -1 


E i*>> 

iG green 

E 1®)’ 

iG brown 


1 


E I®)' 

iGwhite 
















Spatial Search by Quantum Walk with Multiple Marked Vertices 


27 



1 JM-t- 





Vm 


(b) 

Fig. 14 Apart from a factor of — 7 , (a) the Hamiltonian for the second case of search on 
the simplex of complete graphs with k = M + 1 marked vertices, and (b) the leading-order 
terms. 


In this subspace, the search Hamiltonian 0 is 


( 


H = —7 


y/M - 1 
y/M - 1 M — 2 + - 


\ 


1 0 

0 0 

■i VM - 1 

7_ 

y/M - 1 M - 2 
0 0 

0 1 

0 0 


0 ' 
0 
0 
0 

y/M - 2 

y/M - 2 


y/M - 2 yJM - 2 M -2 ) 


This Hamiltonian can be visualized as shown in Figure & The leading-order Hamil¬ 
tonian excludes edges that scale less than y/~M, and it can be visualized as shown in 
Figure E3>- The two eigenstates of that we want to be degenerate are 


y/M + V 8 + M 
with corresponding eigenvalue 


l e ) ' 


y/M + y/8+ M 


I/> + 1 9) 


and 


Eu = -1 (m + y/M(M + 8 )) 

v = — M + 4 — a/m) |a) + | 6 ) 

with corresponding eigenvalue 

E v = i ^—2 - M 7 - y/M y/M + 47 ^ . 

These are degenerate when 7 takes its critical value of 

7 c — ,— . , - , v ~ 1 + — + 0(1/M 2 ). 

\/M (VM + 8 - y/M + 4) M 

The perturbation which restores terms of constant weight, causes certain linear com¬ 

binations a u \u) T cty |u) to be eigenstates of H-|- H^. Doing the perturbative calculation 
to find the coefficients, the perturbed eigenstates for large N are 


(l*> + \v)) , 


E = —3 — 


E = -3- 


1 

y/M 

1 

y/M 
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Since | u) ~ | g) and |u) ~ |6) for large TV, the system evolves from |s) ~ | g) to |6), which is 
marked, in time tt/AE: 

7 T\/ M 
U = ——• 

2 

We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7c = 1 + 3 /M—a straightforward calculation shows that it must be within o(l/M 3 / 2 ). 
This 7 c and runtime are in agreement with Table [2] 


B.3 Two Marked Per Complete Graph, Case 1 


As shown in Figure [T^l, the system evolves in a 2 -dimensional subspace spanned by 

1 


|a) = 

|ft> = 


vmri j 


. = y 1 i). 

V(M + l)(M- 2 ) . e + e 


In this subspace, the search Hamiltonian 0 is 

H= ( 2+1 72(M^2) 

\y / 2 (M - 2 ) M -2 

We can find the eigenvectors and eigenvalues of this directly without perturbation theory. 
They are 

1 + 47 - M7 — vT+^t^tmtT^m^t^, . ... 

= - 2727M^2 7 -' a> + |6> 

with corresponding eigenvalue 

Eq = - ^—1 — M7 + \/l + 87 — 2M7 + M 2 7 2 ^ 


and 


1 + 47 — M7 + y/l + 87 — 2M7 + M 2 7 2 

*01 = - , - \a ) + 6) 

2V2VM - 27 


with corresponding eigenvalue 

E\ = - ^—1 — M7 — \/l + 87 — 2M7 + M 2 7 2 ^ . 
When 7 takes its critical value of 


7c= m’ 


these become for large TV 


IV’o) = + !“>)’ ®o = -i- 


|V> 1> = ^ ( |6 > -|a», Ei = —1 + • 
So the system evolves from |s) « |6) to |a) in time it/AE: 

t\/ M 


t* = 


2 V 2 ' 


An explicit calculation as in m shows that 7 must be chosen within o(l/M 3 / 2 ) of its 
critical value 7 C = 1 /M for this evolution to occur asymptotically. 

This 7 c and runtime are in agreement with Table | 2 | 
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Fig. 15 Apart from a factor of — 7 , (a) the Hamiltonian for the second case of search on the 
simplex of complete graphs with k = 2 (M + 1 ) marked vertices, and (b) the leading-order 
terms. 


B.4 Two Marked Per Complete Graph, Case 2 

As shown in Figure [TJj, the system evolves in a 3-dimensional subspace spanned by 

i a ) = , = y i») 

VmTT)^' 

1 6> = 1 = y |i) 

V(M-4)(Af+l) 

| C > = 1 y i») 

V2( M + !) ie ^w 

In this subspace, the search Hamiltonian § is 

( l+i y/2 (M - 4) 3 \ 


H = - 7 


y/2(M - 4) M - 4 y/2(M - 4) 

V 3 y/2 (M - 4) 1 ) 


This Hamiltonian can be visualized as shown in Figure ET The leading-order Hamil¬ 
tonian excludes edges that scale less than M, and it can be visualized as shown in 

Figure [If- Clearly, the eigenstates of this are |a), 16), and |c) with corresponding eigen¬ 
values — 1, — 7 M, and 0. Since |s) ~ | 6 ), we choose 7 so that \b) and |a) are degenerate, 


7c = 


M 


The perturbation H^\ which restores edges of weight <9(\/M), causes certain linear com¬ 
binations a a \o) + 0 : 5 1b) to be eigenstates of H The coefficients can be found in 
the usual way, and they yield perturbed eigenstates 


IV’o) =9' + 1“))’ 


1 


E 0 ~ —l — \ — 


IV'i) = ~ !“))> 


E\ = —1 - 

So the system evolves from |s) ~ | 6 ) to | a) in time 7r/AE : 


t* = 


We can again use the method of m to find how precisely 7 must be chosen to its critical 
value 7 C = 1/M—a straightforward calculation shows that it must be within o(l/M 3 / 2 ). 
This 7 C and runtime are in agreement with Table [ 2 ] 



















